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Abstract 

We point out that the results of many baryogenesis scenarios operating at or below the TeV 
scale are rather sensitive to the rate of anomalous fermion number violation across the elec- 
troweak crossover. Assuming the validity of the Standard Model of electroweak interactions, 
and making use of previous theoretical work at small Higgs masses, we estimate this rate for 
experimentally allowed values of the Higgs mass (m# = 100. ..300 GeV). We also elaborate 
on how the rate makes its appearance in (leptogenesis based) baryogenesis computations. 
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1. Introduction 



The scenario of thermal leptogenesis relies on anomalous baryon + lepton number vio- 
lation [2j, which is very rapid at temperatures above the electroweak scale jS], to convert 
the original lepton asymmetry into an observable baryon asymmetry. Usually the tempera- 
ture range where the lepton asymmetry generating source terms are active, is much above 
the electroweak scale. In this case the anomalous processes have ample time to operate, 
and their precise rate is not important. In fact, the conversion factors are simple analytic 
functions EI > for which various limiting values were derived already long ago jH [3] . 

However, baryon asymmetry generation may also be a low temperature phenomenon, in 
which CP-breaking source terms are active down to the electroweak scale; for recent examples, 
see Refs. (S]-ES]- I* 1 this case the temperature dependence of the anomalous rate does play 
an important role. This is even more so for the large (Standard Model like) Higgs masses that 
are currently allowed by experiment the electroweak symmetry gets "broken" through 
an analytic crossover rather than a sharp phase transition |18| I19j , whereby the anomalous 
rate also decreases only gradually. 

To allow for a precise study of generic scenarios of this type, it is the purpose of this 
note to collect together all the relevant rate equations, such that systematic errors from this 
part of the computation can be brought under reasonable control. We reiterate the baryon 
and lepton violation rate equations in Sec. [2J estimate the anomalous "sphaleron" rate as a 
function of the Higgs mass and temperature in Sec. |3J and summarise in Sec. HJ 



2. Baryon and lepton number violation rates 

To zeroth order in neutrino Yukawa couplings, the Standard Model allows to define three 
global conserved charges: 

X t = — - U , (2.1) 

n G 

where B is the baryon number, Lj the lepton number of the i th generation, and n G denotes the 
number of generations. Given some values of Xi, a system in full thermodynamic equilibrium 
at a temperature T and with a Higgs expectation value v m i n (suitably renormalised and in, 
say, the Landau gauge), contains then the baryon and lepton numbers [H] 

B = B cq = Xy-7p-J y^ X i > L i = L i,cq = — — - Xi , (2.2) 
i=l 

4[5 + 12n G + 4n G + (9 + 6n G )a " 
65 + 136n G + 44n G + (117 + 72n G )x 2 

These relations hold up to corrections of order 0((Xi/VT s ) 2 ) from the expansion in small 
chemical potentials, 0((hv m i n /TrT) 2 ) from the high-temperature expansion, as well as 0(h 2 ) 
from the weak-coupling expansion, where h is a generic coupling constant. 
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If we deviate slightly from thermodynamic equilibrium, the baryon and lepton numbers 
evolve with time. A non-trivial derivation |3J yields the equations [1J H3 HO] 

B(t) 



B{t) 
p(x) 



-n G p 



nun \ r ( jiff(r) r^^\ 



T 3 



B 



eqj > 



3[65 + 136n G + 44n G + (117 + 72n G )i 



(2.4) 
(2.5) 



2n G [30 + 62n G + 20n G + (54 + 33n G )a; 2 ] 

In the literature the factor n G p(vmin/T) is often replaced with the constant 13n G /4, which 
indeed is numerically an excellent approximation. The term T^isiT) is called the Chern- 
Simons diffusion rate, or (twice) the sphaleron rate, and is defined by 

(Q 2 (t))T 



r d iff(T) 



lim 

V,t-*oo 



Vt 



(2.6) 



where Q(t) = Ldt' J v d 3 x.' q(x') = N cs (t) — N cs (0) is the topological charge, and N cs (t) is 
the Chern-Simons number. The expectation value in Eq. (|2.6[) is to be evaluated in a theory 
without fermions 4 . Corrections to Eq. (|2.5[) are of the same type as those to Eq. (|2.3[) . 

For practical purposes, it is useful to eliminate the conserved charges Xi from the equations, 
and write just a coupled system for B(t), Li{t). Defining 

rdiff(T) _ X{Vmin/T) 



7 = n G p 



1 



X 



T 3 



>1 



1 - X{Vmin/T) ' 



(2.7) 



and introducing sources fi(t) for the lepton numbers, we can convert Eqs. (|2.2j) . (|2.4jl to 



1Q 



B(t) = - 7 (i) B(t)+ V (t)Y, L i(t) 



Li(t) 



i=l 

=i 



(2. 



(2.J 



These equations can easily be integrated, if we know the temperature dependence of Vmm/T 
and the time evolution of T. The solution is particularly simple if we make use of the fact 
that r\ is, to a reasonable approximation, a constant, rj(t) ~ 0.52 ± 0.03. In this case linear 
combinations of Eqs. (|2.8f) . (|2.9|) yield independent first order equations for B(t) — L(t) and 
B(t) + r]L(t), where L(t) = E"=l L i(*)- Denoting w(i';t) = exp[-(l + ??) f*,dt" j(t")] and 
/(*) — Sr=i /*(*)) ^ e solution reads 

>■* 



B(t) 



1 + r? 



B{t Q )+r]L{t ) u{t ;t)+r] B(t )-L(t ) 



-r, J dt'f(t') l-u(t';t) 
'to 



(2.10) 



A further simplification follows by noting that uj(t'; t) varies very rapidly with the time t' 
around a certain t' ~ t*, from zero at t 1 < to unity at t' > t*, while fit') is a slowly varying 
function of time. Assuming furthermore that B(to) = L(to) = 0, we obtain 



Bit) 



-v 



1 + V 



dt'fit') 



to 



X I ' 'dt'fit') 
I to 



(2.11) 
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where the "decoupling time" can be defined as U = t + Jjdi' [1 - u(t'; t)]. Thus, if /(f) / 
around the time i*, the baryon asymmetry generated depends sensitively on i*, and it is 
important to know the function u)(t';t), determined by j(t"), quite precisely. 

The equations that we have written were formally derived in Minkowski spacetime. They 
are easily generalised to an expanding background, however: their form remains invariant if 
we simply replace the total (comoving) baryon and lepton numbers B, by number densities 
over the entropy density s(T): B — > n B = B /[a 3 s(T)], L — > ul = L/[a 3 s(T)], where a 3 is 
a comoving volume element. Furthermore, it is often convenient to replace time derivatives 
with temperature derivatives via 

d _ y/m Vgr) d 

dt mpi d[lns(T)]/dTdT ' { ' ' 

where e(T) is the energy density; we assumed the Universe to be flat (k = 0); and we ignored 
the cosmological constant. Both s(T) = p'(T) and e(T) = Ts(T) — p(T) follow from the 
thermodynamic pressure p(T) which is known to high accuracy |2J, but can in practice be 
reasonably well approximated with the ideal gas formula p(T) « g^T 4 / '90, with ~ 106.75. 

In many baryogenesis scenarios, the source terms fi(t) in Eq. (|2.9jl are approximated by 
Boltzmann-type equations for the various left-handed and right-handed neutrino number 
densities. Collecting the number densities to the matrices n^,n^, respectively, with the 
normalization Tr [n^] = til, a concrete realization of Eqs. (|2.8|) . I|2.9|) could then read 

n B (t) = -7(t){n B (t) + 7?(t)Tr[n L (t)]} , (2.13) 

n L (t) = ~ — {n B (t) + ri(t)Tr[n L (t)]}l + r L [n R ,n L ,t] , (2.14) 
n G I- J 

n R (t) = T R [n R ,n L ,t] , (2.15) 
with functionals TliFr that need to be determined for the specific model in question. 



3. Chern-Simons diffusion rate 

An essential role in the rate equations (|2.13j) - (|2.15[) is played by the function j(t) whose time 
dependence is, via Eq. I|2.7jl . dominantly determined by rdiff(T), defined in Eq. (|2.6[) . We 
now collect together the current knowledge concerning rdifr(T) in the Standard Model. 

At high temperatures (in the "symmetric phase") the Chern-Simons diffusion rate is purely 
non-perturbative, and needs to be evaluated numerically. So-called classical real-time simu- 
lations HU produce T m (T) = (25.4 ± 2.0)a^T 4 [23], where the number 25.4 is in fact the 
value of a function containing terms like la(l/a w ) |24j . at the physical a w . 

At lower temperatures, the rate is traditionally written in the form [22] 
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Here g is the SU(2) gauge coupling, a w = g 2 /Air; u>- might generically be called the dynamical 
prefactor, and is related to the absolute value of the negative eigenvalue of the fluctuation 
operator around the sphaleron solution; A/t r (AAV) ro t are normalisation factors related to the 
zero-modes of the fluctuation operator; k contains the contributions of the positive modes; 
and -Esph is the energy of the saddle-point configuration (the sphaleron) |26j . 

Most of the factors appearing in Eq. (|3.1|) have been evaluated long ago. In particular, -E S ph 
can be found in Ref. [2Z] for the bosonic sector of the SU(2)xU(l) Standard Model, while 
fermionic effects were clarified in Ref. [2HJ. The zero- mode factors and (the naive version of) 
oj- were evaluated in Refs. |291 I30| . while k was determined numerically in Refs. |311 132 j . 

Unfortunately, it as not a priori clear how accurate the corresponding results are. Indeed, 
Eq. (|3.1|) has an inherently 1-loop structure, but it is known from studies of the electroweak 
phase transition that 2-loop effects, parametrically suppressed only by the infrared-sensitive 
expansion parameter O (hT / irv m i n ) , are large in practice [HE] EI]- Moreover, the naive defini- 
tion of lo- through the negative eigenvalue does not appear to be correct |SS] • 

A reliable determination of T^ff can again be obtained by numerical methods, employing 
real-time classical simulations. Of course classical simulations are not exact either, but they 
do contain the correct infrared physics, and should thus only suffer from infrared-safe errors 
of the type mentioned below Eq. (|2.3|) . Thus, classical simulations allow in principle to 
incorporate the dominant higher order effects, as well as a correct treatment of 

In the "broken symmetry phase" , large-scale classical simulations have been carried out in 
Ref. PS]. Unfortunately, they only extend up to Higgs masses around m# = 50 GeV, and 
were only carried out for certain temperatures (there are some results also at larger Higgs 
masses but with less systematics |37j). While we have not carried out any new simulations, 
we do make use of the observation [HE] that the discrepancy between the numerical results, 
and a certain analytical recipe, of the type reiterated below, appears to be independent of 
the Higgs mass. We thus extend the analytical recipe to large Higgs masses, and add to these 
results a (small) constant correction factor, extracted from Ref. |3H| . In practice, the steps 
are as follows: 

(i) We employ the (resummed) 2-loop finite-temperature effective potential V(v ) in Landau 
gauge, as it is specified in Ref. [HE]- Effects of the hypercharge group U(l) need to be taken 
into account only at 1-loop level, as demonstrated in Ref. ^2] ■ The potential is parametrised 
by the zero-temperature physical quantities mw, mz, w-top> Tiff, ot s (mz), Gf\ their values 
(apart from mu) are taken from Ref. |17j . 

We remark that although this is formally a higher order effect, the effective potential does 
depend on the scale parameter Jl of the MS scheme. One may thus consider various choices of 
/2. We follow a strategy similar to Ref. [HU and write V(v) — V(0) = f^dv' dV(v')/dv'\^ = ^ v iy 
where the scale is chosen as p>(v) = A \J 3A e ffW 2 , where A e ff is the scalar coupling of the 
dimensionally reduced theory [HE] and A is a constant. We consider A = 1.0 as the "reference 
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value", while variations in the range A = 0.25. ..4.0 indicate the magnitude of uncertainties. 

(ii) To avoid threshold singularities at small v related to the Higgs and Goldstone masses, 
we replace the exact 2-loop potential by a polynomial fit around the broken minimum: 

Re[v{v) T ; vm =J2b n (v- <uo n + o((v - v min f) , (3.2) 

n=2 

where v = v/T. We carry out the fit in the range v = (0...1.5) v m i n . Only values v < v m \ a are 
needed for the sphaleron solution, but including some larger values allows for a better fit of 
the curvature around the minimum. We have considered other fit forms as well and find that 
the errors introduced through the fitting are insignificant compared with other error sources. 

(iii) We compute the sphaleron energy E sp \jT for this potential. We assume that the 
use of the 2-loop potential rather than the tree-level potential takes care of the factor k in 
Eq. (|3.1|) . which we thus set to unity. At 1-loop level this can to some extent be demon- 
strated explicitly J32j, but what is more important for us is that any possible errors from 
this approximation are compensated for by step (v) below. The effect of the U(l) group is 
treated perturbatively [2E1> which is an excellent approximation [2Zj- We use an effective 
finite-temperature Weinberg-angle tan 2 (6\v)eS ra 0.315 



(iv) We determine the zero-mode factors jV tr , (AAV) ro t and the dynamical factor as 
described in Ref. (HOI) except that every appearance of the tree-level \(h 2 — l) 2 /4g 2 is replaced 
by the 2-loop potential V (hv m i n ) / g 2 s v^ ia . We also determine the effective gauge coupling g c g 
of the dimensionally reduced theory [SB] , and use g e s instead of g in Eq. (|3.1|) . The effect of 
the zero-mode factors and w_ is to effectively decrease E sp ^/T by about 15%, or by 3. ..10 in 
absolute units. 



(v) Finally we add a correction from Ref. [HE]) which we assume to be a constant: 

-(3.6 ±0.6)1 . (3.3) 



r (Mi) = r (i)-(iv) 

1 diff — 1 diflf ex P 



This correction is in most cases subleading compared with those in step (iv), and goes in the 
opposite direction. It may be noted that there is some latitude with respect to which gauge 
is used for the evaluation of the prefactors appearing in Eq. (I3.1|) [23120]) but since is 
gauge-independent, the non-perturbative correction factor compensates for this as well. 

(vi) Finally, since we rely on an extrapolation of the non-perturbative correction factor to 
larger Higgs masses, we assign a generous overall uncertainty to T^m, in the range 
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r diff (T) 



T 4 



2.0 . (3.4) 



This amounts to roughly three times the error in Eq. (|3.3|) . We stress that even though the 
Higgs masses leading to Eq. (|3.3j) are much smaller than we consider, the values of v m \ n /T 
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Figure 1: The temperatures for which specific values ofv m i n /T (in Landau gauge) are reached, 
as a function of the Higgs mass tuh- For v m i n /T = 1.0 we also show the effects of the 



variations m 



top 



174.3 ± 5.1 GeV (dashed lines) and A = 0.25. ..4.0 (dotted lines). 



are similar, and thus the bulk of the effect in Eq. ()3.3[) should still remain intact, at least in 
the physically most plausible range 100 GeV < m# < 200 GeV. 

In Fig. ^ we show the location of the minimum of the 2-loop effective potential. We 
only consider values for which the infrared sensitive expansion parameter hT /7rv m i n remains 
reasonably small. For higher temperatures, the corresponding rate T^ig extrapolates smoothly 
to the symmetric phase value [37] , like standard thermodynamic observables ^2 EH] • 

The rates T^m are displayed in Fig- El with assumed uncertainties of the order in Eq. (|3.4[) . 
For practical applications, we note that in the range 100 GeV < m# < 200 GeV and for T such 
that -ln[r difr (T)/T 4 ] ps 30. ..50, the results can within our uncertainties be approximated by 



In 



r diff (T) 



2^4 



i+j<2 
i,3>0 



- 150 GeVY 
10 GeV ) 



- 150 GeVY 
10 GeV J 



(3.5) 



with the coefficients 



coo 

C20 



39.6 , cio = 3.52 , c i = -7.09 , 
-0.376 , en = 0.421 , c 02 = 0.170 



(3.6) 



Given rdiff(T)/T 4 , we can finally estimate the decoupling time and/or the corresponding 
decoupling temperature T*, needed in Eq. ()2.11|) . In the limit that Tdisi^/T 4 changes very 
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Figure 2: ln[r diff(T) /T ] as a function of the Higgs mass and temperature. The overall error 
is estimated in Eq. jcl.4\ )- The dotted horizontal line indicates the value which all curves 
approach at large T. The values in the range 100 GeV < mu < 200 GeV can be roughly 
approximated by Eq. h'J.5\) . Note that the rate falls off more slowly at large Higgs masses. 



rapidly with T, the solution is given by the equation Uq pTdis{T*)/T* = H(T*), where H(T) 
is the Hubble rate defined through H 2 (T) = 8vre(T)/3m^. Writing 

+ A{T -T,) + 0{{T -T,f) , (3.7) 

corrections to this leading order approximation are of relative order 0{\/ AT*), which accord- 
ing to Eqs. ()3.6j) is in the one percent range, and thus subdominant compared with other 
error sources. The leading order solution is shown in Fig. |31 

Comparing Fig. 01 with Fig. Q it is seen that T* corresponds to values v m \ n /T = 1.0. ..1.2. 
At the same time, the rate of change of TdirT is less abrupt (^4 is smaller) at large Higgs 
masses, and a sudden decoupling is a less precise approximation. This can be seen in Fig. 01 
where the full function 1 — uj(t';t) appearing in Eq. (|2.10fl is plotted. 



In 



diff 



CO 



T 1 



In 



4. Summary and conclusions 

The main contents of this note are the baryon and lepton number rate equations shown 
in Eqs. (|2,7[) - (|2,9[) . as well as the "sphaleron rate" rdiff(T)/T that enters these equations, 
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Figure 3: The solid line indicates the decoupling temperature as defined in the text (assum- 
ing a constant g* ~ 106. 75), with an error band following from changing T ^{T*) /T^ within 
the range of Eq. The dashed lines show the corresponding anomalous rate. 

shown in Fig. [2 and in Eq. Q3.5[) . With this knowledge, and given that the factors %, p, 
r\ are to a fairly good approximation constants, the equations can be integrated in closed 
form, leading to Eq. (|2.10|) . An even simpler estimate for the baryon number generated in a 
given scenario can be obtained from Eq. Q2.11|) . where i* corresponds to the temperature T* 
shown in Fig. |H1 On the other hand, the most precise results can be obtained by integrating 
Eqs. (|2.7j) - (|2.9j) numerically down to temperatures shown in Fig. El All of these equations 
are model-independent in form; the specific model enters through the source terms f{. 

The biggest uncertainties of our estimates for Fdifj(T)/T 4 originate from the fact that 
systematic numerical studies have only been carried out at fairly small Higgs masses |36l I37| . 
If a Standard Model like Higgs particle is found at the LHC, there is certainly a strong 
motivation for repeating the numerical studies at the physical value of the Higgs mass, in 
order to remove the corresponding error source (Eq. (j3.4j) ) from our estimates. 
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250.0 



T' / GeV 

Figure 4: The function 1 — uj(t';t) appearing in Eq. h2.1U\) . as a function of the temperature 
T' corresponding to the time t' (the final moment t is fixed to the point where T = 100 GeV). 
We indicate temperatures instead of times, because this significantly reduces the dependence 
on the constant ~ 106.75, which has non-negligible radiative corrections This figure 

can be used to gauge the accuracy of the sudden decoupling approximation shown in Fig. [3J 
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